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§ ; ABSTRACT 
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(J \ Surprisingly large column densities of H3 have been detected using infrared 

O I absorption spectroscopy in seven diffuse cloud sightlines (Cygnus 0B2 12, Cygnus 

O ; 0B2 5, HD 183143, HD 20041, WR 104, WR 118, and WR 121), demonstrating 

^ I that H;^ is ubiquitous in the diffuse interstellar medium. Using the standard 

^ I model of diffuse cloud chemistry, our Hj]' column densities imply unreasonably 

Tij- ' long path lengths (~ 1 kpc) and low densities (~ 3 cm~^). Complimentary 

^ ' millimeter-wave, infrared, and visible observations of related species suggest that 

O ' the chemical model is incorrect and that the number density of Hj]" must be in- 

I creased by one to two orders of magnitude. Possible solutions include a reduced 

O ■ electron fraction, an enhanced rate of H2 ionization, and/or a smaller value of 

1- ^ • the dissociative recombination rate constant than implied by laboratory ex- 

I ' periments. 

c/3 . Subject headings: infrared: ISM: lines and bands — ISM: clouds — ISM: 
molecules — molecular processes — cosmic rays 

X: : 

I ^Department of Chemistry, Department of Astronomy & Astrophysics, and the Enrico Fermi Institute, 

■ - - ' University of Chicago, Chicago, IL 60637 

^Present address: Department of Astronomy, University of California at Berkeley, 601 Campbell Hall, 
Berkeley, CA 94720; bjmccall@astro.berkeley.edu 

■^National Optical Astronomy Observatories, Tucson, AZ 85726 

^Gemini Observatory, 670 North A'ohoku Place, Hilo, HI 96720 

^Joint Astronomy Centre, Hilo, HI 96720 

^Department of Astronomy, University of Texas at Austin, Austin, TX 78712-1083 

^Faculty of Science, Okayama University, 3-1-1, Tsushima-naka, Okayama 700-8530, Japan 

^Nobeyama Radio Observatory, Minamimaki, Minamisaku, Nagano, 384-1305, Japan 

^Department of Physics, University of Texas at El Paso, El Paso, TX 79968; and McDonald Observatory, 
University of Texas at Austin, Austin, TX 78712 



- 2 - 



1. Introduction 

The molecular ion has long been considered to play an important role in the chem- 
istry of dense molecular clouds, as it initiates the network of ion-neutral reactions (Herbst 
and Klemperer 1973; Watson 1973) that is responsible for the wealth of molecules observed 
by infrared and radio astronomers. However, Hg has not been considered to be an impor- 
tant species in diffuse clouds, because it is thought to be destroyed rapidly by dissociative 
recombination with electrons, which are abundant in diffuse clouds. For the purposes of 
the present study, we do not make a distinction between diffuse and translucent clouds, but 
simply use "diffuse" to indicate sightlines where CO is not the dominant reservoir of carbon 
(in contrast to dense clouds). 

The detection of in the diffuse interstellar medium toward the Galactic Center 
(Geballe et al. 1999) and toward the visible star Cygnus 0B2 12 (McCall et al. 1998a) 
with similar column densities to those of dense clouds (McCall et al. 1999) was therefore 
quite surprising. These observational results imply either a very long pathlength (hundreds 
of parsecs) of absorbing material, or a serious problem with the standard model of diffuse 
cloud chemistry. Various attempts have been made to explain the abundance of H;^ toward 
Cygnus 0B2 12 while preserving the standard model of the chemistry, either by adjusting 
all parameters to optimize (Cecchi-PesteUini and Dalgarno 2000) or by postulating an 
additional ionization source for H2 molecules peculiar to the neighborhood of the Cygnus 
0B2 association (Black 2000). 

We have conducted a small survey of twelve diffuse cloud sources with the aim of de- 
termining whether the sightline toward Cygnus 0B2 12 is unique, or whether high column 
densities of are common in diffuse clouds. As a result of this survey, we have detected 

in seven diffuse cloud sources (Cygnus 0B2 12, Cygnus 0B2 5, HD 183143, HD 20041, 
WR 104, WR 118, and WR 121), and obtained upper limits toward five other sources (HD 
194279, HD 168607, P Cygni, Ori, and C Oph). For Cygnus 0B2 12 and 5 and HD 
181343, we have also obtained infrared and radio spectra of CO, as well as high-resolution 
visible spectra of relevant atoms and molecules. We have also obtained infrared CO spectra 
towards WR 104 and WR 108. 

These results demonstrate the ubiquity of H^ in diffuse clouds, and suggest that there is 
a "global" problem with the current models. In this respect, H^ represents a major problem 
in diffuse cloud chemistry, reminiscent of the Diffuse Interstellar Bands and CH"*". 
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2. Observations and Data Reduction 

The stellar parameters of our targets are listed in Table 1, and a summary of the obser- 
vations is provided in Table 2. The line of B.^ at 3.715 fim was observed using the 
Phoenix spectrometer on the Mayall 4-meter telescope at Kitt Peak National Observatory 
(KPNO). The observations were performed in three runs (July 1998, June 2000, and 
March 2001). The CGS4 facihty spectrometer on the United Kingdom Infrared Telescope 
(UKIRT) was used in July 2000 and May 2001 to study the i?(l,0)-i?(l,l)" doublet of H+ 
near 3.668 //m. The three observed rotation- vibration transitions belong to the funda- 
mental band of the z/2 degenerate bending mode, and arise from the lowest allowed rotational 
states {J,K) = (1,1) and (1,0). The old notation [e.g., -R(l,l)'^] has been superseded by a 
more flexible notation [e.g., i?(l,l)"] which uses u and / in place of + and — [for a complete 
description of the new standard notation for Hg , see McCall (2000) or McCall and Oka 
(2000)]. 

Infrared spectra of portions of the fundamental (f=l-0) and overtone (f=2-0) infrared 
bands of CO were also obtained. The fundamental band was studied toward Cygnus 0B2 
12 and Cygnus 0B2 5 with Phoenix in July 1999 and June 2000, and toward HD 183143, 
WR 104, and WR 118 with CGS4 in May 2001. The overtone band of CO was searched for 
with Phoenix toward HD 183143 in June 1997 and toward Cygnus 0B2 12 in July 1999. 

All the infrared data were reduced using the procedure outlined in McCall et al. (1999), 
which is described in more detail in McCall (2001). Object spectra were ratioed by (scaled) 
standard spectra in order to remove the many strong atmospheric absorption lines in these 
regions. The wavelength calibration was achieved using telluric absorption lines, and is 
estimated to be accurate to roughly 2 km/s. The resolving power of the Phoenix spectrometer 
was roughly 40,000 until 1998, and approximately 60,000 during 1999-2001. The CGS4 
spectrometer was used with the long camera, yielding a resolving power of 40,000. 

High resolution visible spectra of atoms and molecules of interest were obtained using 
the Coude spectrometer on the Smith 2.7-m telescope at McDonald Observatory. A total 

of three echelle settings were used. The flrst covered the K I line at 7699 A and the CN 
A-X v=2-0 i?2(0) line at 7875 A. The second covered the C2 A-X v=2-0 band near 8760 A. 
These two echelle settings were used with the F1E2 conflguration and the 0.2" slit, which 
yielded a resolving power of about 200,000. The third echelle setting covered Ca I at 4227 
A, Ca II at 3934 A, CH at 4300 A, CH+ at 4232 A, and CN B-X v=0-0 near 3874 A. This 
setting was used with the FlEl conflguration and the 0.6" slit, to yield a resolving power 
of approximately 120,000. The McDonald data were reduced using standard IRAF routines 
(CCDPROC and doecslit). It was not generally necessary to ratio the visible spectra by 
standard stars, due to the weakness and paucity of atmospheric absorption lines in the regions 
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of interest. Wavelength calibration was achieved using comparison spectra of a Th-Ar lamp, 
and is estimated to be good to about 0.01 A in all cases. 

Rotational spectra of CO in the Cygnus 0B2 association were obtained at the Nobeyama 
Radio Observatory (NRO) for J— 1-0, and the Caltech Submillimeter Observatory (CSO) 
for J— 2-1 and 3-2. The CSO spectra were obtained using position switching with an "off" 
position of 30' south. This position has significant ^^CO emission, but no detectable ^"^CO 
emission compared with a more distant off position, which is important as the ^^CO data is 
used in the interpretation of the infrared CO lines. The NRO spectra were obtained using 
position switching with an off position of (/ = 81,6 = 3), which was found to be more free 
of CO emission. CO spectra toward HD 183143 were also obtained at the James Clerk 
Maxwell Telescope (JCMT) for J=2-l, and at NRO for J=l-0, using frequency switching. 
In all cases, data were reduced using the standard procedures. 



3. Results 
3.1. Spectra 

In our 1998 observing run at Kitt Peak, we detected for the first time toward Cygnus 
0B2 5, which hes approximately 6' (3 pc at 1.7 kpc) away from Cygnus 0B2 12. Although 
the H.^ R{l,iy line was unfortunately Doppler shifted under the telluric CH4 feature near 
37152 A, we are confident that the detection is secure, because the other telluric CH4 line 
(near 37158 A) is removed very effectively by ratioing with a standard star, and because we 
observed a similar feature in Cygnus 0B2 12 at the same velocity with an equivalent 
width consistent with our earlier observations (McCall et al. 1998a). Our Kitt Peak 
observations are summarized in Figure 1. 

Our 1998 run also yielded detections of Hj]" toward the Wolf-Rayet stars WR 121, WR 
104, and WR 118. The source WR 118 is interesting in that it shows two H3 components, 
which have similar (though not the same) spacing as the two telluric CH4 lines. Consequently, 
the scaling and ratioing by the standard star was performed by choosing parameters which 
minimized the residuals of CH4 lines elsewhere in the spectrum (not shown in the figure) . It 
is worth noting that the inferred velocities of the peaks (~-|-5 and +48 km/s) would occur 
at heliocentric distances of ~0.5 and 3 kpc, respectively, assuming a fiat Galactic rotation 
curve with = 220 km/s, Ro = 8 kpc, and / = 21°. 80. These distances do not exceed the 
estimated distance to the source (m — M = 14, or ~ 6.3 kpc) obtained from the observed 
V magnitude of 22, the visual extinction Ay = 12.8 from van der Hucht et al. (1996), and 
the absolute magnitude My = -4.8 for WC9 stars (Smith, Shara, and Moffat 1990). 
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Finally, our 1998 Kitt Peak run yielded a non-detection toward HD 194279, toward 
which CH is observed with a heliocentric velocity of -11.8 km/s (G. Galazutdinov, private 
communication). This source, along with the Wolf-Rayet stars, was chosen for study because 
the 3.4 fim aliphatic carbon feature (a signature of diffuse clouds) had been detected by 
Pendleton et al. (1994). 

Our June 2000 run at Kitt Peak yielded the first detection of toward a more tra- 
ditional diffuse cloud source, HD 183143. Two components of Hj^ were seen, in agreement 
with the visible spectra of other molecules (see below), although the blue component was 
somewhat affected by the telluric CH4 line. During this run, we also obtained non-dctcctions 
toward P Cygni [which has an LSR velocity of CH"*" of -9 km/s (Hobbs 1973)] and toward ( 
Oph [which has an LSR velocity of CO of -0.79 km/s (Lambert, Sheffer, and Crane 1990)]. 
Our March 2001 run at Kitt Peak yielded a detection of toward the least reddened star 
to date, HD 20041. We also obtained an upper limit of Hj absorption towards O^i. 

Our run at UKIRT in July 2000 (see Figure 2) confirmed the detections of Hj toward 
Cygnus 0B2 5 and HD 183143. In the case of HD 183143, the two-component velocity struc- 
ture is very clear. For both the Cygnus 0B2 5 and HD 183143 observations, instrumental 
artifacts were present in the "negative" spectrum — these features changed in intensity (and 
sometimes disappeared) depending on exactly which rows of the array were extracted, while 
the H^ features were consistently present. We also obtained a non-detection toward HD 
168607 [which has a heliocentric CH velocity of 22 km/s (Gredel 1999)]. 

In May 2001 at UKIRT, we confirmed the detections of Hg" toward the Wolf-Rayet 
stars WR 104 and WR 108 by using the i?(l,l)"-R(l,0) doublet. The separation of the 
two velocity components toward WR 118 is comparable to the spacing between the two Hjj" 
transitions, so that the i?(l,0) of the blue component and the i?(l,l)" of the red component 
are overlapped. 

A table of the observed Hj]" line parameters is given in Table 3. The observed equivalent 
widths have been converted to column densities under the assumption that the lines are 
optically thin, using the transition dipole moments calculated by J. K. G. Watson, as listed 
in McCaU et al. (1998a). 

3.2. CO towards the Cygnus OB2 Association 

A brief look at the high resolution infrared spectrum of CO toward Cygnus 0B2 12 
and Cygnus 0B2 5 was obtained during test and engineering time with Phoenix at Kitt 
Peak in July 1999. FoUowup observations with longer integration time in June 2000 yielded 
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higher signal-to-noise spectra, displayed in Figure 3. The observations of Cygnus 0B2 12 
clearly show the two closely-spaced (~ 5 km/s) velocity components observed in the mm- 
wave spectrum from JCMT (Geballe et al. 1999). The infrared lines suggest a CO excitation 
temperature of ~10 K toward Cygnus 0B2 12, but this should not be taken as an indication 
of the kinetic temperature of the gas due to the effect of spontaneous emission (Geballe et 
al. 1999). The infrared CO line parameters are listed in Table 4. We have adopted the 
transition dipole moments of Hure and Roueff (1996). 

The signal-to-noise ratio of the Cygnus 0B2 12 spectrum is clearly not high enough to 
rely on the measured equivalent widths of the individual components, particularly because of 
the contamination of the telluric lines. However, we can estimate the column density of CO 
by assuming that the true equivalent widths of the two components have the same ratio as 
the integrated areas of the components in the ^^CO emission spectrum, and by adopting the 
6-values of the ^^CO in estimating the saturation corrections. This analysis was performed 
with the high-quality ^^CO spectrum obtained at CSO (shown in Figure 4; the mm-wave 
line parameters are listed in Table 5). We assume that the 7 km/s peak carries 0.39 of 
the equivalent width (and has 6=0.75 km/s), and that the 12 km/s peak carries 0.61 of the 
equivalent width (and has 6=0.88 km/s). With these assumptions we estimate a total CO 
column density of ~ 1.4 x 10^^ cm~^ in front of Cygnus 0B2 12. This estimate is considerably 
higher than the estimate of 3 x 10^^ cm~^ given in Geballe et al. (1999), but still much lower 
than the total column density of carbon (~ 2.5 x 10^^ cm"^) inferred from the color excess. 

We have also obtained an upper limit on the CO column density toward Cygnus 0B2 
12, based on our failure to detect the f =2-0 overtone with Phoenix. The upper limit on the 
v—2-0 features yield an upper limit {3a) of N(CO) < 3 x IQ-'^''' cm~^. This suggests that at 
most about 10% of the carbon along the line of sight to Cygnus 0B2 12 can be in the form 
of CO, and that this sightline does not consist of dense clouds. 

Using the ^^CO data from CSO for J=2-l and from NRO for J=l-0 (see Figure 5), we 
performed an analysis using an LVG code to calculate the statistical equilibrium, assuming 
a temperature of about 30 K (consistent with the C2 analysis of §3.5 and the H3 analysis of 
§4.2). If the density {n — 200 cm~^) derived from the analysis of C2 toward Cygnus 0B2 12 is 
assumed, the ^^CO column density, iV(13) is determined uniquely by matching the observed 
line emission. The results are iV(13) = 1.8 x 10^^ cm~^ for the 7 km/s component and 
iV(13) = 2.7 X 10^^ cm-2 for the 12 km/s component. The sum of these, A^(13) = 4.5 x 10^^ 
cm~^, would imply a total column density of CO of 2.3-4.5 x 10^^ cm^^ for isotope ratios 
of 50-100; this column density is about twice that inferred from the infrared CO absorption 
lines. If, on the other hand, we enforce agreement with the column density from the infrared 
absorption, and assume an isotope ratio of 100, the ^^CO data indicate somewhat higher 
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densities: n — 500 cm~^ for the 7 km/s component, and n — 900 cm"^ for the 12 km/s 
component. While these densities are somewhat higher than inferred from C2, it would not 
be surprising to have density variations, with the ^^CO data probing the denser gas, and the 
C2 probing somewhat less dense gas. 

The hne of sight toward Cygnus 0B2 5 appears to have considerably less CO than that 
toward Cygnus 0B2 12, based both on the marginal detections in the infrared spectrum 
(Figure 3) as well as the CSO ^^CO spectra (Figure 4). Therefore we expect that this 
sightline is also dominated by diffuse cloud material, although because of the low signal-to- 
noise of the infrared spectrum it is difficult to estimate the CO column density. 

In order to probe the spatial extent of the CO gas in the Cygnus 0B2 association, 
we have obtained (position-switched) ^^CO J=l-0 spectra at four other locations in the 
association, labelled A through D in Table 2 and Figure 6, which also shows the spectra 
toward Cygnus 0B2 12 and 5. 

3.3. CO towards HD 183143 

To check for the presence of CO in this sightline, we obtained rotational spectra of the 
J=2-l line (at JCMT) and the J=l-0 line (at NRO). The JCMT results (frequency-switched 
at 8 MHz and 16 MHz) are shown in Figure 7. While the J=2-l spectrum of Cygnus 0B2 
12 showed nearly 2 K of emission at the H3 velocity, the spectrum of HD 183143 shows 
no emission at the correct velocities in excess of ~ 0.1 K. The NRO results are shown in 
Figure 8, and provide an upper limit of ~ 0.2 K for the J=l-0 emission at the observed Hj 
velocities. The closest reasonable velocity component is at +25 km/s, and is only about 0.5 
K. In contrast, Cygnus 0B2 12 shows 3 K of emission. 

During our May 2001 UKIRT run, we obtained an absorption spectrum of the i?-branch 
of the CO I'^l-O fundamental band. This spectrum (see Figure 9 and Table 4) shows weak 
absorptions near vi^sR ~ 25 km/s, in agreement with the radio spectra. Because of the fairly 
small (geocentric) Doppler shift (~ 11 km/s) at the time of the observation, it is difficult to 
accurately measure the column density, but it is clear that the total CO column density is 
less than 10^^ cm~^. A previously obtained Phoenix spectrum of the v=2-0 overtone band 
(from June 1997) provides a less stringent upper limit on the total CO column density of 
- 10^^ cm-2. 

Both the radio and infrared results demonstrate that there is very little (if any) CO 
associated with the H3 toward HD 183143. 
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3.4. CO towards WR 104 and WR 118 

At UKIRT in May 2001, we obtained spectra of the CO fundamental band toward the 
Wolf-Rayet stars WR 104 and WR 118. The spectra are shown in Figure 9 and the hne 
parameters are hsted in Table 4. WR 104 shows a single velocity component centered at 
I'LSR ~ 21 km/s — surprisingly, this is different from the velocity of the H3 (~ 10 km/s)! 
WR 118 has a very complicated velocity structure, which we have fit with four Gaussian 
components. The lowest velocity component at vi^sr ~ 10 km/s may be associated with the 
H3 component at ulsr ~ 5 km/s, and the blended components near i>lsr ~ 40-55 km/s are 
consistent with the second H3 component at vlsr ~ 48 km/s. 

Without any radio spectra of these sources, it is difficult to attempt saturation correc- 
tions, so we have simply listed in Table 4 the column densities in the hmit that the hues are 
optically thin. In this limit, WR 104 has N(CO) ~ 9 x 10^^ cm-^ and WR 118 has N(CO) 
~ 4 X 10^^ cm-2 (in the J=0-3 levels). 

3.5. C2 Spectra 

The C2 A-X v—2-0 band was clearly detected toward Cygnus 0B2 12 at McDonald, as 
shown in Figure 10. Each line appears as a clear doublet with a separation comparable to 
that of the CO, except for the Q{8) line which is blended with P(4). The equivalent widths 
were measured separately for the two velocity components, except for the Q{8) + -P(4) blend, 
for which only the total equivalent width could be measured. The contribution of P(4) was 
estimated from the strength of Q(4), which permitted an estimate of the Q{8) equivalent 
width for the sum of the two velocity components. The (5(10) line was contaminated by an 
atmospheric H2O absorption line, and the spectra were ratioed by that of a Cep in order 
to remove the contribution of the telluric line. The derived C2 line parameters (for each 
component separately, as well as for the total profile) for Cygnus 0B2 12 are hsted in Table 
6. Column densities have been determined assuming the lines are optically thin, using the 
oscillator strength /20 = 1.67 x 10"^ of van Dishoeck (1983). 

Using the method of van Dishoeck and Black (1982), and assuming the scahng factor 
for the radiation field I — 1 and the C2-H2 coUisional cross-section ao — 2 x 10~^^ cm^, we 
can estimate the number density n of collision partners, as well as the kinetic temperature 
T. This analysis is performed by calculating the rotational distribution for a grid of points 
on the (n.T) plane, then finding the best fit to the column densities N{J) derived from the 
observed spectrum. The advantage of this method is that it yields the range of values of 
(n,T) that are consistent with the observations, given the uncertainties in the measurements. 
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A web-based "C2 calculator" is available at littp://dib.ucliicago.edu/c2 

For the 7 km/s component, we obtain a best fit value of (n,T) = (220 cm~^, 40 K), and 
for the 12 km/s component we obtain (210 cm~^, 30 K). Using the total column densities 
rather than those of the individual components yields (200 cm~^, 30 K). These values are 
fairly well constrained by the observational data — it is exceedingly unlikely (i.e. it requires 
more than a 3a error in at least one of the measurements) that the density is outside the 
range 150-600 cm~^ or that the temperature is outside the range 25-55 K. 

The spectrum of Cygnus 0B2 5 also shows a hint of the doublet structure, but the 
signal-to-noise of the spectrum is not sufficient to reliably measure the equivalent widths of 
the individual components. Therefore we have measured only the total equivalent widths, 
which are reported in Table 7. An analysis of the rotational excitation (while very uncertain) 
suggests that the density exceeds 800 cm~^ and that the temperature lies in the range 60-90 
K. According to Gredel and Miinch (1994), chemical models predict a sudden decrease in C2 
abundance at densities above 10^-^ cm~^. Taken with our (uncertain) results, this suggests 
that the density in the region where C2 is found is in the range 800-3000 cm~^. 

In agreement with the results of Grcdcl (1999), wc saw absolutely no trace of C2 towards 
HD 183143. Assuming the Q{2) line would have a width of ~ 5 km/s, the 3a limit on the 
equivalent width is about 1.7 mA, which corresponds to a column density of N( J=2) < 3x 10^^ 
cm~^ (comparable to Gredel's upper limit), which is more than ten times less than that of 
Cygnus 0B2 12! 

3.6. Other Visible Spectra 

The results of our high resolution visible spectroscopy are shown in Figure 11 (HD 
183143), Figure 12 (Cygnus 0B2 5), and Figure 13 (Cygnus 0B2 12). These figures are 
plotted in velocity (with respect to the local standard of rest), and also display the infrared 
measurements of and CO. 

CH appears as a doublet in Cygnus 0B2 5, with velocities consistent with the infrared 
CO measurements. CH also appears as a doublet in HD 183143, with remarkably similar 
velocity structure to the Hj. The blue wavelength (4300 A) of the CH transition made it 
inaccessible in the case of Cygnus OB 12 (which is very heavily reddened), but it may be 
possible to detect the line with a larger telescope such as the Hobby-Eberly Telescope. The 
line parameters of CH (as well as CH"*" and CN) are listed in Table 8. The CH column 
densities have been derived using the curve of growth of van Dishoeck and Black (1989), 
assuming 6=1 km/s. 
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It is interesting to compare the observed CH column densities with those which would 
be predicted from the empirical relation between H2 and CH (Magnani and Onello 1995): 

Ar(H2) = 2.1 X lO^iV(CH) + 2.2 x 10=^° cm'^ 

For HD 183143, E(B-V)=1.2, so Nh ~ 7x 10^^ [assuming the standard gas-to-dust conversion 
factor (Bohlin, Savage, and Drake 1978)], and (if / = 2N(H2)/[N(H)+2N(H2)] = 2/3) N(H2) 
~ 2.3 X 10^^ cm^^. From the empirical relation, one would then predict N(CH) ~ 1.0 x 10^^, 
roughly twice the observed vahie. For Cygnus 0B2 5, with E(B-V) = 2.1, we estimate 
N(H2) ~ 4 X 10^^ and therefore N(CH) ~ 1.8 x 10^^ cm~^, over three times the observed 
value. Assuming our choice of 6-value is not too large, this suggests that these two sources 
may be somewhat different chemically from the usual diffuse cloud sources used to develop 
the empirical relation. However, this discrepancy may not indicate a drastic departure in 
chemical conditions, as other sources with lower color excess have been observed with similar 
departures from the empirical relations (Danks, Federman, and Lambert 1984). The CH 
discrepancy can be reduced by adopting a lower value of /, but this would decrease ^(Hg ), 
exacerbating the problems discussed in §4.2. 

Another interesting molecule is CN: it is observed (marginally) as a doublet in both 
Cygnus 0B2 12 (A-X) and in Cygnus 0B2 5 (A-X and B-X). However, only one component 
(at vlsr ~ 25 km/s) is observed toward HD 183143 — whereas Hj^, CH, and CH+ all 
show two velocity components! The R(l) line is also detected toward HD 183143, and the 
population of the J = and J = 1 levels are consistent with the temperature of the cosmic 
microwave background. 

Finally, it is worth noting that, at least in the case of Cygnus 0B2 12 (Figure 13), the 
velocity profile of H^ appears more similar to that of K I than to C2, CN, or CO. This could 
be because these three molecules are concentrated in the denser regions of the sightline, 
whereas H^ (and K i) may exist over a larger path length. 

4. Discussion 

4.1. Model of Chemistry in Diffuse Clouds 

The formation of Hj, in dense or diffuse clouds, begins with the (cosmic ray induced) 
ionization of H2 to form Hj. The generally assumed rate of ionization is C ~ 3 x 10"^^ s~^, so 
that the average H2 gets ionized roughly once every 10^ years. Once ionized, the H^ quickly 
reacts with another H2 to form H^ and an H atom. This ion-neutral reaction proceeds with 
the Langevin rate constant of 2 x 10~^ cm^ s~^, so that in a medium of H2 density 100 cm~^. 
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the average must wait about 2 months to react to form H3 . Clearly the initial ionization 
is the rate-limiting step, so we can say that H3 is formed at a rate of C n(H2). 

The destruction of is very different in diffuse clouds, compared to dense clouds. 
In dense clouds, the dominant destruction path is an ion-neutral reaction with CO, with a 
rate constant of ~ 2 x 10~^ cm^ s~^. In diffuse clouds, electrons are very abundant (due 
to photoionization of C), so dissociative recombination dominates, with a rate constant of 
ke = 4.6 X 10"^ cm^ s'^ (Sundstrom et al. 1994). 

Although the electrons produced by photoionization are probably formed at a high 
temperature (photons are available at energies up to 13.6 eV, and the ionization potential of 
C is only 11.3 eV), they will be quickly thermalized by collisions with H2. Such collisions will 
occur with a Langevin rate constant (which scales as where ^ is the reduced mass) 

of ~ 8 X 10^^ cm^ s~^; in a medium with n(H2) = 100 cm~^, an electron experiences such a 
collision about once per day. The lifetime of the electron is limited by radiative recombination 
with C"*" ions, which occurs with a rate constant of ~ 2 x 10~^^ cm^ s~^ (Aldrovandi and 
Pequignot 1973) at the low temperatures (~30 K) of diffuse clouds. If n(C"'") = 10~^ cm~^, 
the average electron lifetime is about 10^ years. Evidently the electrons will be very quickly 
thermalized to the kinetic temperatures of the clouds. The electron recombination rate 
constant is then /cg ~ 5 x 10~^ cm^ s~^, and the destruction rate of H3 can be given as kg 
n{Rt) n{e). 

To calculate the number density of (McCall et al. 1998a; Geballe et al. 1999), we 
make the steady-state approximation that the rates of formation and destruction are equal. 
Thus, 

Cn(H2) = ken(RtMe) 



which can be rearranged to give 



^ ke n{e) 



The chemistry of H3 is unique in the sense that the number density of H3 is not 
dependent upon the number density of the cloud, but only on the ratio of the number 

densities of molecular hydrogen and electrons (as well as the ratio of two constants C/^e)- 
We can further simplify the above equation by considering /, the fraction of protons in H2, 
writing n(H2) = (//2)[2n(H2) +n(H)] = (//2)n(SH). Since most electrons are formed from 
the ionization of carbon, we can also write n(e) = (1 — «)n(SC), where a is the fraction 
of carbon in un-ionized (neutral or molecular) form. Finally, the gas-phase ratio of carbon 
atoms to hydrogen atoms is usually defined as zq = n(SC)/n(SH). Putting this all together. 



- 12 - 



we find 

^(H+) l ^— — (1) 

^ ^' ke 2 1-a zc ^ ^ 

Adopting tlie values C = 3 x 10^^'' s^^, = 5 x 10"'' cm^ s~-^, / = 2/3, a ~ 0, and 
zc = 1.4 X 10"^ (Cardelli et al. 1996; Sofia et al. 1997), we obtain 

n(H+) = 1.4 X 10~^ cm-^ (2) 

The validity of the steady-state approximation can be checked by comparing the timescale 
needed to reach steady-state with other relevant timescales. The timescale needed to achieve 
steady-state is approximately the steady state number density of Hj]" divided by the forma- 
tion rate, or n{}i^) /[(^n(H2)], which for n(H2) = 100 cm~^ is about one year — clearly much 
shorter than other relevant timescales! 



4.2. Inferred Cloud Parameters 

One cloud parameter can be determined from the H3 observations independent of the 
chemical model — the kinetic temperature of the gas. This is possible because the ortho 
(J=l, K=0) and para (J=l, K=l) levels of are efficiently thermalized together through 
proton hop and proton exchange reactions with H2 (McCall et al. 1998b). These reactions 
occur with the Langcvin rate constant of 2x 10~^ cm^ s~^, so that the average experiences 
such a reaction about every two months (assuming n{B.2) = 100 cm^^). The lifetime of the 
average H3 can be estimated from the dissociative recombination rate ke and the number 
density of electrons n{e) ~ 10~^ cm~^ to be about 4.5 years, considerably longer than the 
thermalization timescale. Therefore, we can use the Boltzmann expression 

^ortho _ 9ortho ^-AE/kT _ 2g-32.87/r 
^para 9para 

to estimate the kinetic temperature from the observed ortho.-para ratio (in this equation, the 
g values are the statistical weights). For Cygnus 0B2 12, we obtain T=27±4 K (McCall et 
al. 1998a); for Cygnus 0B2 5, 47±13 K; for HD 183143, 31±3 K; for WR 104, 38±10 K; 
and for the -|-48 km/s component of WR 118, 40±3 K. [We cannot estimate the ortho.-para 
ratio for the low velocity component toward WR 118 due to the overlapping of the lines.] 

Given our calculated number density of from the chemical model, we can now 
estimate the path length of absorption, using the relation L — N (H'^) / nCR^) . Once the 
path length has been calculated, we can then estimate the average number density of collision 
partners along the path length as (n) — [A'"(H2) -|- A'"(H)]/L, where A'"(H2) and A'"(H) can be 
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estimated from the color excess and an assumed value of / (2/3). The results of this analysis 
are given in Table 9. 

In most cases, the derived pathlengths are a substantial fraction of the estimated dis- 
tance to the star, which seems difficult to accept. In addition, the derived average number 
densities (for the sources in which Hg was detected) are in the range 1-5 cm~^, which seems 
unreasonably low. These densities are nearly two orders lower than the densities typically 
derived from the rotational excitation of C2. In addition, these densities are so low that a 
substantial fraction of H2 should be photodissociated, meaning that our value of / = 2/3 
should be lowered, which would in turn make the number density even smaller, the path 
lengths even longer, and the average density even lower. 

It seems clear that there is a serious problem with the model, and that n{}i^) must be 
larger (probably by at least one order of magnitude) than we have calculated. In the next 
subsection, we explore possible solutions to this problem. 



contains five parameters: the H2 ionization rate {(), the dissociative recombination rate 
[kg), the fraction of protons in molecular form (/), the fraction of carbon atoms that are 
not ionized (a), and the gas phase carbon fraction (^c)- We now consider each one of these 
parameters in more detail. 

The fraction of protons in molecular form (/) has been assumed to have the value 
2/3, which is the largest value found in the studies of molecular and atomic hydrogen by 
Copernicus, lUE, and FUSE. The maximum value of / is 1, so even making this adjustment 
(as is done in the model of Cecchi-Pestellini and Dalgarno (2000)) will only increase ^^(H^) 
by a factor of 1.5 — hardly enough to fix the problem. In addition, increasing the value of / 
would increase the discrepancy between inferred and measured A'^(CH), as discussed in §3.6. 

The gas phase carbon fraction (^c) has been taken to be 1.4 x 10^'^ based on ultraviolet 
observations of classical diffuse clouds (Cardelli et al. 1996; Sofia et al. 1997). Lowering this 
value (for instance, by assuming that a larger fraction of the carbon is depleted onto grains) 
would reduce the number density of electrons and therefore increase the number density 
of H^. However, there is no astronomical evidence to support the magnitude of depletion 
necessary to mitigate the Hg problem. 



4.3. 



Possible Solutions 



The equation for the number density of H3" 




(3) 
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The fraction of carbon that is not ionized (a) has been assumed to be near zero. It is 
conceivable that a substantial fraction of the carbon could be in the form of neutral C atoms 
(rather than C^), and there is no observational evidence (in these particular sightlines) that 
requires most of the carbon to be ionized. However, most chemical models of translucent 
clouds (see, for example. Figure 9 of van Dishoeck and Black (1988)) suggest that when the 
optical depth has increased to the point where n(C)=n(C"'"), over 10% of the carbon atoms 
are already in the form of CO. Taken along with our observations of a low column density of 
CO in these sighthnes, this suggests that a must be less than 0.5, and thus could contribute 
at most a factor of 2 towards increasing the Hg number density. 

However, other models (e.g. Fcderman and Huntress (1989)) suggest that the abundance 
of C"*" in an Ay — 4 cloud may be 15 times lower than in an Ay — 1 cloud. If we adopt 
this model, then it is possible that a <; 0.9, which would increase ^(Hg ) by one order of 
magnitude. Given the differences between the models, we hesitate to rely on them too 
much, in the absence of observational evidence. The best determination of the value of a 
would come from an observation of the C"*" transitions at 1334 and 2325 A using the Hubble 
Space Telescope. 

Along the same lines, one might speculate about the possibility of an "electron sink" in 
these sightlines, so that C could still be ionized but the free electron abundance would be low. 
The best candidate mechanism for removing electrons from the gas is probably attachment 
to grains or large molecules (e.g., Lcpp and Dalgarno (1988)). However, because the number 
density of large molecules or grains is orders of magnitude lower than the number density 
of electrons (since these molecules or grains form from elements with cosmic abundance less 
than or equal to that of carbon, and each molecule or grain contains many atoms), this 
process could not effect the removal of a significant fraction of the electrons from the gas 
phase. 

The H2 ionization rate {Q has been assumed to be 3 x 10~^^ s~\ but this value is not 
terribly well constrained. Cecchi-Pestellini and Dalgarno (2000) assume 6 x 10~^^ in order 
to increase the Hg number density, and values as high as ~ 2 x 10~^^ have been derived from 
analysis of the chemistry leading to OH in diffuse clouds (van Dishoeck and Black 1986). 
While the flux of high energy (> 100 MeV) cosmic rays can be constrained by observations in 
the interplanetary medium, the flux of lower energy cosmic rays is essentially unconstrained, 
due to the influence of the Sun. If the cosmic ray spectrum is assumed not to roll off too 
rapidly below 100 McV, it is conceivable that ( might be as high as 10~^^ s~^ (Hayakawa, 
Nishimura, and Takayanagi 1961). A large flux of low energy cosmic rays might increase 
the HJ number density in diffuse clouds while not seriously affecting the situation in dense 
clouds, into which the low energy cosmic rays could not penetrate. 
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Other sources of H2 ionization have also been suggested. Black (2000) has suggested that 
X-rays from the luminous stars in the Cygnus 0B2 association could increase the effective ( 
and therefore n(Hj]"). While this suggestion might solve the problem for the sightlines toward 
Cygnus 0B2, our observations of high iV(Hj]") toward several other sources imply that a more 
general solution is needed. Black (2000) has also suggested that ultraviolet photoionization 
of H2 might contribute. However, photoionization of H2 requires photons with energies above 
15.4 eV, which will ionize H atoms. While Black suggests that such photons could escape 
the H II region near the Cygnus 0B2 giants, it seems unlikely that they could penetrate the 
boundaries of diffuse clouds, where H atoms are abundant. Furthermore, our new detections 
of again suggest that the sightlines toward the Cygnus 0B2 association are not unique. 

The last parameter in the equation for the H3 number density is the dissociative re- 
combination rate (ke). In this work, we have adopted the value of ke — 4.6 x 10~^ j'-o.es 
cm^ s~^ derived from storage ring experiments (Sundstrom et al. 1994). However, the value 
of this constant has been the matter of great controversy over the past two decades — other 
experimental techniques currently yield a value of that is about one order of magnitude 
lower (for a review of the field, see Larsson (2000)). It has been suggested that the dis- 
crepancy might be due to stray fields present in the storage ring experiments, and that the 
recombination rate under interstellar conditions could be quite low. To make matters worse, 
attempts to theoretically calculate the recombination rate yield rates more than an order of 
magnitude smaller than the smallest values obtained in experiments (Orel, Schneider, and 
Suzor-Weiner 2000), but it should be noted that the theory of this recombination process is 
not yet mature. Given the present uncertainty in the true value of k^, it is possible that the 
problem of diffuse cloud Hj may be solved on this front. 



4.4. Observational Tests 

The resolution of the mystery of diffuse cloud Hj will most likely come from further 
observational work. First we must determine with certainty whether the unexpectedly high 
number density of H3 is a "local" problem (i.e., due to special conditions in these particular 
sightlines, or this class of sightlines), or if it is a "global" problem. The present observational 
situation is summarized in Figure 14, which plots the observed Hjj' column density (or the 
upper limit, denoted by a bar to zero) versus the color excess E(B-V). In addition to the 
present observations, the detection of H^ toward the Galactic Center source GC IRS 3 
(Geballe et al. 1999) has been added (the H3 column density plotted is that of the narrow 
component attributed to diffuse clouds by (Geballe et al. 1999), and the adopted color excess 



is that attributed to diffuse clouds by (Whittet et al. 1997)). Keeping in mind that 



7V(H+) ~ n(H+) L oc n{R+) 



E{B-V) 



(4) 



the plot is not inconsistent with a constant ?T.(Hjj') and small variations (of a factor of a few) 
in the average density (nn)- Therefore at this stage there is no strong evidence that the 
"local" parameters (/, a, and zc) are any different in the detection sources than they are in 
the non-detections such as C Oph. 

However, this inference needs to be tested directly by observations. This can be achieved 
by detecting along a line of sight with lower E(B-V) which can be studied by FUSE and 
HST. Once the column densities of H3 , H I, H2, C I, C ll, CO, and CH are obtained for a 
single sightline, the influence of the "local" parameters will be directly determined. 

Assuming that these parameters are found to have their canonical values, we will be left 
with the ratio (/kg. Absent a speedy resolution from the experimenters or theorists in the 
field of dissociative recombination, further insight into these two constants can perhaps best 
be obtained through an observational survey of in more heavily reddened lines of sight. 
The rate of destruction is controlled by n(e)/nH (which depends on the optical depth in 
the ultraviolet), whereas the rate of Hg formation is controlled by C (which, with a given 
incident cosmic ray spectrum, depends on the stopping power of the cloud as a function of 
cosmic ray energy). Therefore, a detailed study of the diffuse-to-dense cloud transition, along 
with chemical models, could help constrain the values of ( and ke- The highly-reddened, 
early-type subset of the Stephenson catalog recently compiled by Rawlings, Adamson, and 
Whittet (2000) may serve as a good starting point for such a study. 



The H3" molecular ion has now been definitely detected in seven diffuse cloud lines of 
sight, suggesting that its unexpectedly high abundance is not due to the peculiarities of a 
particular region (the Cygnus 0B2 association), but rather a general feature of the diffuse 
interstellar medium. is observed in clouds with and without CO, C2, and CN, confirming 
that the chemistry that leads to is completely decoupled from that which is responsible 
for these heavier diatomics. The most hkely explanation for the high abundance is a 
larger than expected C/K ratio — due either to a larger flux of low-energy cosmic rays or 
to a lower value of ke in interstellar conditions. The possibility of a lower than expected 
electron density has not been ruled out, and should be directly tested by finding Hg in a 
less reddened source. 



5. 



Conclusions 
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Fig. 1. — Spectra of the B.^ line observed with Phoenix at Kitt Peak. Each object 

spectrum has been ratioed by that of a standard star to remove atmospheric absorption Hnes. 
The vertical arrows denote the observed (thick arrows) or expected (thin arrows) positions 
of Hg depending on each source's Doppler shift at the time of observation. 
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Fig. 2.— Spectra of the H+ i?(l,0)-i?(l,l)" doublet of H+ observed with CGS4 at UKIRT. 
Each object spectrum has been ratioed by that of a standard star to remove atmospheric 
absorption hnes. The vertical arrows denote the observed or expected (for HD 168607) 
position of the Hg doublet. The crosses label instrumental artifacts. 
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Fig. 3. — Infrared v—1-0 CO spectra obtained with Phoenix at Kitt Peak. Object spectra 
have been ratioed by standard star spectra to remove the atmospheric CO hnes. The hori- 
zontal segments of the spectra indicate regions where the strong telluric absorptions could 
not be ratioed out. 
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Fig. 4. — CO spectra obtained at CSO. 
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Fig. 6. — CO J=l-0 "map" of the Cygnus 0B2 association obtained at NRO. Spectra have 
(K) as vertical axis, vlsr (km/s) as horizontal axis. All coordinates are J2000. 
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Fig. 7.— ^^CO J=2-l frequency switched spectra toward HD 183143, obtained at JCMT. 
The velocities of the H3 hnes are marked with vertical dashed hnes. 
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Fig. 8. — ^^CO J=l-0 frequency switched spectra toward HD 183143, obtained at Nobeyama 
(folding has not been performed). The velocities of the lines are marked with vertical 
dashed lines. 
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Fig. 9.— ^^CO ^;=l-0 fundamental band observations of HD 183143, WR 104, and WR 118, 
obtained at UKIRT. Object spectra have been ratioed by standard star spectra to remove 
the atmospheric CO hnes. 
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Fig. 10. — Spectra of the A-X v—2-0 band of C2 obtained at McDonald. 
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Fig. 11. — Summary of spectra of HD 183143, in velocity space. 
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Fig. 12. — Summary of spectra of Cyg 0B2 5, in velocity space. The structure near 30 and 
60 km/s in the lower trace is due to atmospheric lines. 
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Fig. 13. — Summary of spectra of Cyg 0B2 12, in velocity space. The structure near 30 and 
60 km/s in the lower trace is due to atmospheric lines. 
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Table 1: Stellar Parameters. 



Name 


V 


E(B-V) 


Spectral Type 


Distance (pc) 


WR 118 


22.0 


4.13^ 


WC9 


6300^ 


Cyg 0B2 12 


11.4 


3.35^ 


B5Iab 


1700^ 


WR 104 


13.5 


2.10^ 


WC9 


1300^ 


Cyg 0B2 5 


9.2 


1.99^ 


07e 


1700<^ 


WR 121 


11.9 


1.68^ 


WC9 


1690"^ 


HD 168607 


8.25 


1.6P 


B9Ia 


1100^ 


HD 183143 


6.9 


1.28^ 


B71a 


looof 


HD 194279 


7.1 


1.22^ 


B1.51a 


iioof 


HD 20041 


5.8 


0.70S 


AOIa 


1400^ 


P Cygni 


4.8 


0.63*^ 


B2pe 


1800^ 


Ori 


4.6 


0.44' 


B2Iae 


1000^ 


COph 


2.6 


0.32^ 


09V 


140j 



"From Pendleton et al. (1994), assuming Ry = Ay / E(B-V) = 3.1. 
''Smith, Shara, and Moffat (1990). 
=From Schultc (1958). 

'^Torres-Dodgen, Carroll, and Tapia (1991). 
'^Prom Snow, York, and Welty (1977). 
■''Our estimate (spectroscopic parallax). 
sRacine (1968). 

''From Lamers, de Groot, and CassatcUa (1983). 
*From the intrinsic color of Wcgner (1994). 
•'From Hipparcos catalog (Ferryman et al. 1997). 
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Table 2: Log of observations. 



Date (UT) 


Observatory 


Species 


Transition 


Object 


Standard 


Time (min) 


Jun 25, 1997 


KPNO 


CO 


v=2-Q 


HD 183143 


OL Lyr 


60 


Jul 3, 1998 


KPNO 


H+ 

3 


R(l,l)' 


Cye OB2 12 


Q: Cyg 


30 


Jul 3, 1998 


KPNO 


H+ 

3 


-R(l,l)' 


Cve OB2 5 


OL Cve 


164 


Jul 4, 1998 


KPNO 


H+ 


R(l iV 


WR 104 


c\ Tivr 


35 


Jul 4, 1998 


KPNO 


H+ 

3 


R(l lY 


WR 121 


o/ TjVT" 


112 


Jul 4, 1998 


KPNO 


H+ 

^^3 


R(l iV 


WR 118 


Ca Over 


46 


Jul 4, 1998 


KPNO 


H+ 


R(l iV 


HD 194279 


>-* 


92 


Jul 4, 1999 


KPNO 


CO 


v=2-0 


CvK OB2 12 


(X Cyg 


4 


Jul 5, 1999 


KPNO 


CO 


v=l-0 


Cyg OB2 12 


a Lyr 


4 


Jul 5, 1999 


KPNO 


CO 


v=l-0 


Cyg OB2 5 


OL Lyr 


8 


Jun 25, 2000 


KPNO 


H+ 

3 


R(l,lY 

J. j-y 


CvK OB2 12 


a Cyg 


60 


Jun 25, 2000 


KPNO 


H+ 

3 


J. oi^ J-, j-y 


HD 183143 


OL Lyr 


84 


Jun 25, 2000 


KPNO 


H+ 

3 


R(l lY 


P Cygni 


CK Lyr 


78 


Jun 26, 2000 


KPNO 


H3+ 


i?(l,l)' 


COph 


a Lyr 


48 


Jun 26, 2000 


KPNO 


CO 


v=l-0 


Cyg OB2 12 


OL Lyr 


36 


Jun 26, 2000 


KPNO 


CO 


v=l-0 


Cyg OB2 5 


Tivr 


90 


Mar 13, 2001 


KPNO 


H+ 
^^3 


R(l lY 


HD 20041 


a CMa 


100 


Mar 13 2001 


KPNO 


H+ 

^^3 


R(l lY 


v2 Ori 


a CMa 


60 


Jul 6, 2000 


UKIRT 


H+ 


R{ifl),R{i,iy 


HD 183143 




29 


Tul fi 2000 


UKIRT 


H+ 

^^3 


R(\ 0) R(l 1)'" 


Cve OB2 5 


Cyg 


43 


Jul 7, 2000 


UKIRT 




i?(l,0),7i(l,l)- 


HD 168607 


TT Sgr 


25 


May 24, 2001 


UKIRT 


H3+ 


_R(l,0),i?(l,l)'^ 


WR 104 


17 Oph 


8 


AAn^f 9/1 9001 


UKIRT 


H+ 
"3 


0~l f}(A 1 


WR lis 

VV IX -LlO 


OL Cyg 


« 
u 


May 28, 2001 


UKIRT 


CO 


V = l-0 


HD 183143 


a Aql 


24 


May 28, 2001 


UKIRT 


CO 


v=l-0 


WR 104 


Sco 


3 


Mav 9nni 


UKIRT 


CO 
f K I 


v=l-0 

2p _ 2 c ^ 


WR IIS 


r? Oph 


5 


Sep 12, 2000 


McDonald 


1 CN 


A V 11—9 


Cyg OB2 12 


a Cyg 


60 


■^pn 1 9 9nnn 




55 




C\r0- OR9 


^ ^y§ 


60 


Spn 19 9000 


MrDnrialH 


55 


„ 


HD 1 S'HI 4*^ 


Cyg 


60 


^pn 1 9 9nnn 


A/TrDnnalH 


Co 


A Y 11—9 


r'vCT OR9 1 9 


a Cep 


60 


^pr> 1 9 9nnn 


ivid-zoiiam 


Co 


A-X 1!— 9-0 


rivo- OR9 ^ 


a ep 


60 


Spn 12 9000 


MrDnnalH 


Co 
02 

• Qa, I 
Ca II 


A-X ?)— 2-0 

2p _ 2^ 


HD 18^14*^ 


a ep 


60 


Spn 14 2000 


MrDnnalH < 


CH 

CH+ 

CN 


A-X ?)— 0-0 
R-X 71-0-0 


. rjvff OB2 "i 


„ 

^ ^yg 


200 


^pn ^A 9nnn 

OtJJJ -L^) ^UUU 


ivici-'oiiaiu 








„ 

^ Cyg 


100 


Tnl 9R 1QQQ 


cso 




J=2-l 


rivo- nR9 1 9 


30' S 


14 


Till 9S 1 QQQ 


cso 


i^CO 


J=2-l 


pvo- nR9 ^1 


30' S 


14 


Till 9S 9Q 1 QQQ 


cso 


13qq 


J=2-l 


r'vo- nR9 1 9 


30' S 


32 


Till 9S 9Q 1 QQQ 


cso 




J=2-l 


Pvcr OR9 ^ 


30' S 


32 


Till 9K 1 QQQ 


cso 


13qq 


J=3-2 


Pvo- OR9 1 9 


30' S 


20 


Tnl 98 1QQQ 


cso 


13 QO 


J=3-2 


Hva^ OR9 ^ 


30' S 


20 


Apr 9, 2000 


NRO 


12 CO 


J=l-0 


Cyg OB2 12 


(i=81,ft=3) 


18 


Apr 9, 2000 


NRO 


"CO 


J=l-0 


Cyg OB2 12 


(;=81,b=3) 


37 


Apr 11, 2000 


NRO 


12CO 


J=l-0 


Cyg OB2 "A" 


(«=81,b=3) 


10 


Apr 11, 2000 


NRO 


12CO 


J=l-0 


Cyg OB2 "B" 


(i=81,b=3) 


10 


Apr 11, 2000 


NRO 


12CO 


J=l-0 


Cyg OB2 "C" 


(i=81,f)=3) 


15 


Apr 11, 2000 


NRO 


12CO 


J=l-0 


Cyg OB2 "D" 


(i=81,6=3) 


10 


Apr 14, 2000 


NRO 


12 CO 


J=l-0 


Cyg OB2 5 


(i=81,6=3) 


13 


Oct 13, 2000 


JCMT 


12 CO 


J=2-l 


HD 183143 


(8 MHz) 


60 


Oct 13, 2000 


JCMT 


12CO 


J=2-l 


HD 183143 


(16 MHz) 


60 


Dec 14, 2000 


NRO 


12CO 


J=l-0 


HD 183143 


(3 MHz) 


12 


Dec 14, 2000 


NRO 


12CO 


J=l-0 


HD 183143 


(6 MHz) 


13 


Dec 15, 2000 


NRO 


12CO 


J=l-0 


HD 183143 


(4 MHz) 


31 


Dec 15, 2000 


NRO 


12CO 


J=l-0 


HD 183143 


(9 MHz) 


45 
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Table 3: Hg Line Parameters. 



Object 


Line 


vlsr 


FWHM 


Wa 


C7(Wa)^ 

(A) 


N(H+) 

/ -1 1 zl 9\ 

(10^* cm ^) 


cr(N)^ 






(km/s) 


(km/s) 


(A) 


(10" cm" 


Cyg 0B2 12 


i?(l,l)' 


7.0 


14.6 


0.044 


0.002 


2.02 


0.09 


Cyg 0B2 5 


R{i,iy 


7.3 


11.5 


0.022 


0.003 


0.99 


0.15 




i?(i,i)" 


5.6 


15.2 


0.035 


0.004 


1.43 


0.16 




i?(l,0) 


5.8 


15.4 


0.048 


0.004 


1.21 


0.10 


WR 121 


i?(l,iy 


9.7 


17.8 


0.024 


0.004 


1.12 


0.20 


WR 104 


i?(l,l)' 


10.9 


15.1 


0.028 


0.002 


1.27 


0.07 




i?(l,l)" 


11.2 


14.6 


0.030 


0.005 


1.23 


0.21 




i?(l,0) 


7.9 


18.5 


0.041 


0.006 


1.05 


0.14 


WR 118 


i?(l,l)' 


4.8 


9.5 


0.019 


0.001 


0.89 


0.06 






47.6 


12.3 


0.041 


0.002 


1.88 


0.07 




i?(l,l)" 


3.7 


17.1 


0.045 


0.003 


1.86 


0.11 






44.3*' 


15.4 


0.026 


0.002 


1.06 


0.10 




i?(l,0) 


8.5'' 


21.7 


0.077 


0.003 


1.95 


0.07 






48.1 


19.6 


0.065 


0.003 


1.64 


0.07 


HD 183143 


i?(l,l)' 


8.3^= 


10.5 


0.020 


0.002 


0.94 


0.08 






23.6 


8.4 


0.024 


0.002 


1.10 


0.07 






7.1 


12.2 


0.012 


0.001 


0.50 


0.05 






23.4 


9.9 


0.017 


0.001 


0.70 


0.05 




R{1,0) 


5.6 


14.0 


0.022 


0.001 


0.54 


0.03 






22.1 


11.0 


0.023 


0.001 


0.59 


0.03 


HD 20041 


R{i,iy 


-0.5 


11.5 


0.038 


0.002 


1.74 


0.09 


HD 194279 


i?(i,iy 




lO'^ 




0.005 




0.2P 


P Cygni 


i?(i,iy 




lO'^ 




0.002 




0.10^ 


COph 


i?(i,iy 




lO'^ 




0.001 




0.06<= 


HD 168607 


^(1,1)" 




10<^ 




0.003 




0.13^ 




i?(l,0) 




lO'^ 




0.003 




0.07^= 


Ori 






lO'^ 




0.003 




0.12<^ 



"Statistical uncertainties (Ict) arc listed. Systematic errors (for instance, due to ratioing of telluric lines) are 
difficult to estimate and arc likely larger. 

''These two features are overlapped, so the fit results are highly uncertain. 
^Affected by telluric line. 

Adopted FWHM for purposes of computing upper limits. 
'^The firm upper limit for N(H3 ) should be considered to be 3a. 
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Table 4: Infrared ^^CO ^=1-0 Line Parameters. 
Object 



Line 


vlsr 


FWHM 


Wa 


ct(Wa)'' 


N(thin)'' 


(7{NY 




(km/s) 


(km/s) 


(A) 


(A) 


(10^'"^ cm"2) 


(10^^ cm 


P2 


5.9 


4.3 


0.199 


0.021 


2.28 


0.24 




11.8 


6.1 


0.300 


0.025 


3.42 


0.29 


PI 


6.0 


4.1 


0.344 


0.024 


4.71 


0.32 




10.3^= 


1.8 


0.127 


0.015 


1.74 


0.21 


RO 


7.5 


6.0 


0.598 


0.029 


2.75 


0.13 




13.1 


3.8 


0.285 


0.023 


1.31 


0.11 


Rl 


6.6 


4.6 


0.437 


0.021 


3.01 


0.15 




12.0 


3.9 


0.355 


0.019 


2.45 


0.13 


R2 


7.6 


5.1 


0.248 


0.010 


1.91 


0.08 




13.1 


4.9 


0.185 


0.010 


1.42 


0.08 




4.9 


4.8 


0.303 


0.041 


2.09 


0.29 




12.2 


4.5 


0.204 


0.041 


1.40 


0.28 


RO 


25.6 


13.0 


0.138 


0.013 


0.64 


0.06 


Rl"^ 


24.4 


16.1 


0.133 


0.010 


0.91 


0.07 


RO 


21.4 


12.9 


0.520 


0.024 


2.39 


0.11 


Rl 


21.2 


12.4 


0.588 


0.019 


4.05 


0.13 


R2 


21.5 


7.2 


0.247 


0.014 


1.89 


0.11 


R3 


20.6 


5.1 


0.143 


0.011 


1.16 


0.09 


RO 


10.5 


18.2 


1.289 


0.020 


5.93 


0.09 




45. O*' 


18.0 


1.147 


0.020 


5.28 


0.09 




59.4« 


15.6 


0.198 


0.018 


0.91 


0.08 




85.0 


14.0 


0.107 


0.017 


0.49 


0.08 


Rl 


9.2 


16.6 


0.962 


0.015 


6.63 


0.10 




41.6^ 


15.1 


1.078 


0.014 


7.43 


0.10 




56.9<^ 


10.5 


0.204 


0.012 


1.41 


0.08 




81.6 


15.1 


0.211 


0.015 


1.45 


0.10 


R2 


9.1 


12.8 


0.248 


0.013 


1.90 


0.10 




41. 7« 


10.6 


0.539 


0.012 


4.13 


0.09 




53.6^= 


13.2 


0.111 


0.013 


0.85 


0.10 




82.9 


11.1 


0.144 


0.012 


1.10 


0.09 




36.4 


19.1 


0.219 


0.010 


1.77 


0.08 




82.3 


12.5 


0.097 


0.008 


0.78 


0.06 



Cyg 0B2 12 



Cyg 0B2 5 
HD 183143 
WR 104 

WR 118 



"Statistical uncertainties {la) are listed. Systematic errors (for instance, due to ratioing of telluric lines) are 
difficult to estimate and may be larger. 

''Listed column densities arc for the lower state of each absorption line in the given velocity component. These 
are direct calculations from the corresponding value of Wa, in the optically thin limit. For a more detailed 
analysis of the CO column density towards Cygnus 0B2 12, see section 3.2. 

■^Bad fit due to interference of telluric line. 

'^Other lines too marginal to fit. 

'^Blend of two components; individual fit parameters uncertain. 
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Table 5: mm-wave CO toward Cygnus 0B2 12. 



Species 


Line 


vlsr 


HWHM 


(peak) 


Area 






(km/s) 


(km/s) 


(K) 


(K km/s) 




J=l-0 


6.93 


1.30 


0.29 


0.52 






12.58 


1.32 


0.41 


0.58 




J=2-l 


6.9 


1.39 


0.21 


0.31 






12.4 


1.35 


0.33 


0.47 




J=l-0 


7.31 


2.52 


1.49 


3.34 






12.72 


1.81 


3.08 


5.98 




J =2-1 


7.48 


1.92 


0.92 


1.72 






12.48 


2.04 


1.74 


3.66 



Table 6. C2 Line Parameters (components labelled as 1 and 2) for Cygnus 0B2 12. 



Line 




N(J)i 




FWHMi 


W| ^ 


N(J)2 




FWHM2 


(mA) 


N(J)*°* 




(mA) 


(IOI2 cm-2) 


(km/s) 


(km/s) 


(mA) 


(IOI2 cm-2) 


(km/s) 


(km/s) 


(IOI2 cm-2) 


R(0) 


4.96 ± 0.60 


4.38 ± 0.53 


6.7 


2.8 


6.87 ± 0.46 


6.06 ± 0.41 


12.6 


1.7 


11.83 ± 0.76 


10.4 ± 0.7 


Q{2) 


9.91 ± 0.83 


17.5 ± 1.5 


6.7 


3.1 


10.67 ± 0.70 


18.8 ± 1.2 


12.4 


2.2 


20.58 ± 1.09 


36.3 ± 1.9 


0(4) 


8.97 ± 0.75 


15.8 ± 1.3 


6.6 


3.2 


9.65 ± 0.59 


17.0 ± 1.0 


12.3 


2.0 


18.62 ± 0.95 


32.8 ± 1.7 


P(2) 


2.05 ± 0.60 




6.8 


2.0 


3.38 ± 0.55 




12.6 


1.7 


5.43 ± 0.81 




Q(6) 


7.40 ± 0.66 


13.0 ± 1.2 


6.6 


3.1 


7.33 ± 0.53 


12.9 ± 0.9 


12.2 


2.0 


14.73 ± 0.85 


25.9 ± 1.5 


/ Q{8) 1 


















19 ± 2 b 




I P{^) J 




















P(4)'= 


















6 ± 0.3 




0(8)'^ 


















13 ± 2 


22.3 ± 5.2 


Q(10) 


3.9 ± 0.6 


6.6 ± 1.2 


(6.7)- 


(3)« 


3.0 ± 0.5 


5.3 ± 1.5 


(12.5)'' 


{2.ir 


6.9 ± 0.8 


11.9 ± 2.0 



''Statistical uncertainties (Icr) are listed unless otherwise noted. Systematic errors are difficult to estimate and may be larger. 

''Estimate of systematic uncertainty. 

'^Estimated from Q(4). 

'^Blend equivalent width minus P(4) estimate. 
"^Values constrained in the fit to the ratioed spectrum. 
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Table 7: C2 Line Parameters for Cygnus 0B2 5. 



Line 




N(J) 




(mA) 


(10^2 cm-2) 


m 


6.6 ± 0.8 


5.8 ± 0.7 


Q(2) 


10.7 ± 1.0 


18.9 ± 1.8 


Q(4) 


13.5 ± 1.0 


23.8 ± 1.8 


Q(6) 


8.0 ± 1.2 


14.1 ± 2.1 


Q{8) \ 


11 ± 3^ 




.^'(4)/ 




P(4)'= 


4.5 ± 0.3 






7 ± 3*^ 


12.3 ± 5.3 




1.5 ± 0.5 


2.6 ± 0.9 



"Statistical uncertainties {la) are listed unless otherwise noted. Systematic errors are difficult to estimate 
and may be larger. 
''Estimate of systematic uncertainty. 
''Estimated from Q(4). 

"^Blend equivalent width minus P(4) estimate. 
'^Marginal feature; values very uncertain. 
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Table 8: Line Parameters for CH, CH+f and CN!^ 



Object 


Species 


Line 


VLSR 

(km/s) 


FWHM 

(km/s) 


Wa 

(mA) 


N = 
(10^2 cm-2) 


Cyg 0B2 5 


CH 


4300 


5.6 


6.0 


26.5 ± 0.7 


44.0 ± 1.2 d 








12.7 


3.4 


7.1 ± 0.5 


9.2 ± 0.6 ^ 


HD 183143 


CH 


4300 


7.7 


7.7 


12.3 ± 0.7 


17.0 ± 1.0 








23.6 


4.9 


18.6 ± 0.6 


27.7 ± 0.9 


Cyg 0B2 5 


CH+ 


4233 


2.6 


5.0 


11.2 ± 0.6 


12.8 ± 0.7 








7.6 


4.1 


7.1 ± 0.6 


8.1 ± 0.7 








13.3 


8.9 


11.8 ± 0.8 


13.5 ± 0.9 


HD 183143 


CH+ 


4233 


9.0 


4.9 


18.3 ± 0.5 


21.0 ± 0.6 








24.8 


4.5 


22.1 ± 0.5 


25.3 ± 0.6 


Cyg 0B2 12 


CN 


A-X R2(0) 


6.8 


1.1 


1.3 ± 0.4 


3.1 ± 1.0 








12.4 


1.7 


4.3 ± 0.6 


10.3 ± 1.4 


Cyg 0B2 5 


CN 


A-X R2(0) 


5.6 


2.0 


1.2 ± 0.3 


2.9 ± 0.7 








13.7 


1.7 


1.4 ± 0.3 


3.4 ± 0.7 






B-X R(0) 


5.6 


3.9 


4.P 


0.9^ 








12.1 


3.9 


2.8'' 


0.6° 


HD 183143 


CN 


B-X R(0) 


25.2 


2.5 


5.6 ± 0.4 


1.3 ± 0.1 






B-X R(l) 


24.1 


1.5 


1.2 ± 0.2 


0.4 ± 0.1 



"For CH+, we adopt the oscillator strength /oo = 5.45 x 10 ^ of Larsson and Sicgbahn (1983). 
''For CN, we adopt the oscillator strengths = 3.38 x 10~^ (Larsson, Siegbahn, and Agren 1983) and 
/2AX = 7 5 X 10-4 (Davis et al. 1986). 

^Statistical uncertainties (Icr) are listed unless otherwise noted. Systematic errors are difficult to estimate 
and may be larger. 

'^Estimated using the CH curve of growth of van Dishoeck and Black (1989), assuming 6=1 km/s. All other 
column densities computed assuming the lines are optically thin. 

^These values are extremely uncertain due to the low flux of Cygnus 0B2 5 at violet wavelengths. 
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Table 9: Cloud Parameters Based on Conventional Chemical Model. 



Object 


iV(H3+) 


L 


d 


E(B-V) 


7V(H2) + iV(H) 






(IQi^ cm-2) 


(pc) 


(pc) 


(mag) 


(10^1 cm-2) 


(cm^'^) 


Cyg 0B2 12 


3.8*^ 


905 


1700'^ 


3.35= 


12.95 


4.8 


Cyg 0B2 5 


2.6<^ 


629 


1700*^ 


1.99° 


7.69 


4.1 


HD 183143 


2.3^ 


552 


1000° 


1.28*' 


4.95 


3.0 


HD 20041 


3.5s 


833 


1400*^ 


0.70"^ 


2.71 


1.1 


WR 121 


2.2s 


524 


1690' 


1.68j 


6.50 


4.1 


WR 104 


2.3<^ 


547 


1300' 


2.10j 


8.12 


4.9 


WR 118 




1548 


6300' 


4.13j 


15.97 


3.4 


HD 194279 


<1.2S 


<286 


llOO'' 


1.22^ 


4.72 


>5.5 


HD 168607 


<0.6 


<143 


1100° 


1.61*^ 


6.23 


>14.5 


P Cygni 


<0.6s 


<143 


1800"^ 


0.63"^ 


2.44 


>5.7 


yd Ori 


<0.7s 


<171 


1000^ 


O.44I 


1.70 


>3.3 


COph 


<0.4s 


<86 


140™ 


0.32f 


1.24 


>4.8 



"Prom McCall et al. (1998a). 

''Torrcs-Dodgcn, Carroll, and Tapia (1991). 
'^From Schultc (1958). 

''The UKIRT observations of i?(l,0) and i?(l,l)" are adopted for these sources. 
^Our estimate (spectroscopic parallax). 
■''Prom Snow, York, and Welty (1977). 

^Total H3 column density estimated to be twice that of •pwra-lS^. 
''Racine (1968). 

*Smith, Shara, and Moffat (1990). 

■JProm Pendleton ct al. (1994), assuming Ry = Ay I E(B-V) = 3.1. 
'^From Lamers, de Groot, and Cassatella (1983). 
'From the intrinsic color of Wegner (1994). 
""Prom Hipparcos catalog (Ferryman et al. 1997). 



